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Abstract 

We compute the hyperfine splitting in a heavy quarkonium composed of different 
flavors in next-to-leading logarithmic approximation using the nonrelativistic renormal- 
ization group. We predict the mass difference of the vector and pseudoscalar charm- 
bottom mesons to be M(B*) - M(B C ) = 50 ± 17 (th) ±\ 5 2 (5a s ) MeV. 

PACS numbers: 12.38.Bx, 14.65.Fy, 14.65.Ha 

1 Introduction 

The recently discovered charm-bottom heavy quarkonium completes the well investigated 
charmonium and bottomonium families and offers a new perspective in the study of the 
nonrelativistic dynamics of the strong interactions. The first experimental observation of 
about twenty events interpreted as the decays of the B c meson by CDF collaboration pQ 
does not match the precision of the spin one charmonium and bottomonium measurements. 
The statistics, however, is expected to increase significantly in future experiments at Tevatron 
and LHC greatly improving the accuracy of the data. Note that only the pseudoscalar (spin 
singlet) state has been observed so far while the vector (spin triplet) meson B* is still to be 
discovered. This distinguishes cb quarkonium from the bb system, where it is the pseudoscalar 
rjb meson, which asks for experimental detection. 

From the theoretical point of view, the charm-bottom mesons are "in between" the 
approximately Coulomb bb mesons and the cc mesons. Therefore, a simultaneous analysis 
of all three quarkonia could shed new light on the balance between the perturbative and 
nonperturbative effects and further check whether a perturbative analysis provides a reliable 
starting point for them. Moreover, since the nonperturbative effects in the cb system are 
suppressed with respect to the cc meson, the former could be a cleaner place to determine 



the charm quark mass (provided the experimental accuracy is good enough). Another point 
to be stressed is that, though the leading order dynamics of the cb state is quite similar to 
the bb and cc one (up to the value of the reduced mass) the higher order relativistic and 
perturbative corrections are different. Thus the comparison of cb and bb (cc) properties could 
help to establish fine details of the nonrelativistic dynamics. 

The spectrum of the charm-bottom quarkonium has been subject of numerous investi- 
gations based on potential models [21 [3] , lattice simulations [I] , and pNRQCD [5] . This last 
analysis computed the ground state energy within a pure perturbative approach. We con- 
sider that this analysis further indicates that a perturbative approach can be a good starting 
point for studing the B c system. 

In the present paper we focus on the hyperfine splitting (HFS) E^ s of the B c , i.e. the 
mass difference between the singlet and triplet spin states M(B*) — M (B c ). The QCD study 
of the heavy quarkonium HFS has a long history [6j [7] . For the same-flavor quarkonium the 
next-to-leading order (NLO) 0(a s ) correction to the ground state HFS can be found in [5] 
in a closed analytical form. The leading order HFS is proportional to the fourth power of 
the strong coupling constant a s (v) and thus the low order calculations suffer from strong 
spurious dependence on the renormalization scale u, which essentially limits the numerical 
accuracy of the approximation. Hence, the proper fixing of the normalization scale becomes 
mandatory for the HFS phenomenology. The dynamics of the nonrelativistic bound state, 
however, is characterized by three well separated scales: the hard scale of the heavy quark 
mass m, the soft scale of the bound state momentum mv, and the ultrasoft scale of the bound 
state energy mv 2 , where v oc a s is the velocity of the heavy quark inside the approximately 
Coulomb bound state. To make the procedure of scale fixing self-consistent one has to resum 
to all orders the large logarithms of the scale ratios. For the same-flavor case this problem has 



been solved in Ref. [H] within the nonrelativistic renormalization group (NRG) approach and 



the next-to-leading logarithmic (NLL) result for HFS has been derived. The renormalization 
group improved result shows better stability with respect to the scale variation. Moreover, 
the use of the NRG significantly improves the agreement with the experimental value of 
HFS in charmonium in comparison to the NLO computation. Below we generalize the 
analysis to the different-flavor quarkonium case and apply the result to predict the splitting 
M{B*)-M{B C ). 

2 Renormalization group running of the spin depen- 
dent potential 

To derive the NRG equations necessary for the NLL analysis of the HFS, we rely on the 
method based on the formulation of the nonrelativistic effective theory [10] known as poten- 
tial NRQCD (pNRQCD) [11]. The method was developed in Ref. [12] where, in particular, 
the leading logarithmic (LL) result for HFS has been obtained (see also Ref. [13]). In pN- 
RQCD the HFS is generated by the spin-flip potential in the effective Hamiltonian, which 
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in momentum space has the form 



8H apin = D$ - — - — S 2 , S = — - , (1) 

■ 3mim 2 2 

where <7\ and <x 2 are the spin operators of the quark and antiquark with masses m\ and 
m 2 , Cp = (N 2 — l)/(2iV c ), and Dg? f g is the Wilson coefficient, which incorporates the effects 
of the modes that have been integrated out. In effective theory calculations such couplings 
become singular as a result of the scale separation. The renormalization of these singularities 
allows one to derive the equations of the NRG, which describe the running of the effective- 
theory couplings, i.e. their dependence on the effective-theory cutoffs. The solution of these 
equations sums up the logarithms of the scale ratios. 

In general, one should consider the soft, potential and ultrasoft running of corre- 
sponding to the ultraviolet divergences of the soft, potential, and ultrasoft regions [14] . We 
denote the corresponding cutoffs clS V g , Up and u us , respectively. v us and v p are correlated as 
was first realized in Ref. [15]. A natural relation between them is v us = i/|/(2m r ), where 
m r = mxm^l (mi + m 2 ) is the reduced mass. The dependence on u s first emerges in the LL 
approximation after integrating out the hard modes. It disappears after subsequent integrat- 
ing out the soft modes giving rise to a dependence on k, the three-dimensional momentum 
transfer between the quark and antiquark. Thus the soft running effectively stops at v s — k. 
The dependence on v p emerges for the first time in the NLL approximation and cancels out 
in the time-independent Schrodinger perturbation theory for heavy quarkonium observables. 

(2) 

Thus, in pNRQCD one considers D s2 ' as a function of k and v p . For the calculation of the 
spectrum it is convenient to expand this /c-dependent potential around k = v s 

Dfl(k, u p ) = D${v» u p ) + In (If) v a -£-D$(u a , v p ) + . . . . (2) 

yl/g J ell's 

The characteristic momentum for the Coulomb system is a s m r and for v s ~ a s m r the average 
of In (k/ u s ) over bound state wave function does not produce a large logarithm while the 
derivative in In v s results in extra factor of a s . Thus, for the calculation of the HFS in NLL 
approximation one can take the first term on the right hand side of Eq. ([2]) in the NLL 
approximation, the second term in the LL approximation and neglect the higher derivative 
terms. 

Once expanded, the potential is a function of v s and v p (we should not forget that there 
is also a dependence on rrii, the masses of the heavy quarks, and Uh, the matching scale of 
the order of the heavy quark masses). Let us start with the discussion of the soft running. 
To the NLL approximation it is determined by the following NRG equation 

Us ~d~v~ D ^' s = a s c F(m 1 )c F (m 2 )-f s , (3) 
where cf is the effective Fermi coupling, 

7T 7T Z 
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is the soft anomalous dimension and a s = a s (u s ) is renormalized in the MS scheme. The 
running of the coefficient cp is known in NLL approximation [16]. It reads 
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where z = (a s (z/ s )/a s (z// l )) 1///3 °, ^ ~ m« is the hard matching scale, c\ = 2(Ca + Cp) and the 
one- and two-loop anomalous dimensions read [16] 
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Here Ca = N c , Tp — 1/2, in is the number of massless quark flavors, and /3j is the (i + l)-loop 
coefficient of the QCD (3 function 



A) = y Ca - ^ T ^ n ' ' A = 
The value of one-loop anomalous dimension 
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can be extracted from the result of Ref . [12] . The result for the two-loop coefficient 
1 
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is new. It was obtained by an explicit calculation of the subleading singularities of the 
two- loop soft diagrams using the approach of [TTJ [HI [19]. In this approach, dimensional 
regularization with D = 4 — 2e is used to handle the divergences, and the formal expressions 
derived from the Feynman rules of the effective theory are understood in the sense of the 
threshold expansion [H]. Thus the practical calculation reduces to the evaluation of the 
coefficients of the quadratic and linear soft poles in e. Our approach possesses two crucial 
virtues: the absence of additional regulator scales and the automatic matching of the con- 
tributions from different scales. For the reduction of the two-loop Feynman integrals to the 
master ones the method of Ref. [2D] was used. 

The solution of Eq. Q can be written as a sum of the LL and NLL contributions. The 
LL result is already known and reads [12] (see also [T3] ) 
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For the NLL term we obtain 
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The potential running starts to contribute in NLL order. To compute it we inspect all opera- 
tors that lead to spin-dependent ultraviolet divergences in the time-independent perturbation 
theory contribution with one and two potential loops [12J, [2T], [22] . They are 

(i) the (D(v 2 ,a s v) operators [5J, 

(ii) the tree 0(v 4 ) operators, some of which can be checked against the QED analysis 



(iii) the one-loop 0(a s v 3 ) operators for which only the Abelian parts are known [T5], while 
the non- Abelian parts are new. 

In the NLL approximation, we need the LL soft and ultrasoft running of the 0(v 2 ) and 0(v 4 ) 
operators, which enter the two-loop time-independent perturbation theory diagrams, and the 
NLL soft and ultrasoft running of the 0(a s v) and 0(a s v 3 ) operators, which contribute at 
one loop. The running of the 0(v 2 ,a s v) operators is already known within pNRQCD [12J. 
The running of the other operators is new. For some of them, it can be obtained using 
the reparameterization invariance [24] . We refrain from writing the corresponding system of 
NRG equations, which is rather lengthy, and only present its solution, which can be cast in 
the form 
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where the coefficients A{ and fi are given in the Appendix. To get this result we rescale the 
ultrasoft cutoff to v us = v^jvh. The difference to the previous definition is beyond the NLL 
accuracy. 

The LL result ffTUl) obeys the tree level matching condition 



LL 



Ot S {Vh) 



(15) 



while Eqs. ( II 1|) and ffl4l) vanish at v — Vh by construction. We then use the known one-loop 
result of the potential [6J to obtain the NLO matching condition at the scale k = v s = u p = v^. 
It reads 
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Note that in the limit mi = m 2 = m q this equation does not reproduce the same-flavor 
equal-mass expression 
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because of the two-gluon annihilation contribution present in the latter case. 
Thus the NLL approximation for the Wilson coefficient is given by the sum 
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where = (z/, z/) and we combine the soft and potential running by setting 

v s — v v — v, which is consistent at the order of interest. From Eqs. ([1]) and (jSJ) we obtain 
the final result for the NLL spin-flip potential 
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3 Hyperfine splitting in NLL approximation 

We are now in the position to derive the NLL result for the HFS. It is obtained by computing 
the corrections to the energy levels with the insertion of the potential ffl9|) in the quantum 
mechanical perturbation theory. The result for principal quantum number n reads 
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where v = 2C F a s m r , E° = -C 2 F a 2 s m r / (2n 2 ) , # ft (z) = d n \nT{z)/dz n , T(z) is the Euler T- 
function, and 7^ = 0.577216 ... is Euler's constant. In Eq. ( 1201) the first order correction to 
the Coulomb wave function at the origin due to one-loop contribution to the static potential 
reads [25] 
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where a% = 31Ca/9 — 20T F rii/9. Furthermore, the second line of Eq. (1201) results from the 
second term in square brackets in Eq. (TT9~j) after average over the Coulomb wave function. 
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By expanding the resummed expression up to O(o%), we get 
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where a s = a s (z/), \I/„(x) = <i n lnT(x)/(ix n , L" s = In (Cpo; s /n) and 2^ = 2m r and v = u/n 
has been chosen. The C(a 2 ln 2 a s ) term is known [T2], [13], while the 0(a 2 lna s ) term is 
new. The equal- mass case expression [9J, relevant for charmonium and bottomonium, can 
be deduced from Eq. (j20|) by replacing 
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and setting mi = m 2 . After including the one-photon annihilation contribution, the Abelian 
part of the equal-mass result reproduces the (D(ma® In a s ) and 0(ma 7 s \n 2 a s ) corrections to 
the positronium HFS (see e.g. [23l [T8] ) . 



4 Numerical estimates and conclusions 

For the numerical estimates, we adopt the strategy of [9] and replace the on-shell mass of 
the charm and bottom quarks by one half of the physical masses of the ground state of 
bottomonium and charmonium [26]. In practice, we take m& = 4.73 GeV and m c = 1.5 
GeV, consistent with the accuracy of our computation. Furthermore, we take a s (Mz) as an 
input and rurQ with four-loop accuracy down to the matching scale Vh to ensure the best 

1 For the running and decoupling of a s we use the program RunDec [27] . 
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precision. Below the matching scale the running of ct s is used according to the logarithmic 
precision of the calculation in order not to include next-to-next-to-leading logarithms in our 
analysis. In Fig. HJ the HFS for the charm-bottom quarkonium ground state is plotted as 
a function of v in the LO, NLO, LL, and NLL approximations for the hard matching scale 
value Vh = 2.05 GeV. As we see, the LL curve shows a weaker scale dependence compared 
to the LO one. The scale dependence of the NLO and NLL expressions is further reduced, 
and, moreover, the NLL approximation remains stable at the physically motivated scale of 
the inverse Bohr radius, Cpa s m r ~ 0.9 GeV, where the fixed-order expansion breaks down. 
At the scale v' 0.81 GeV, which is close to the inverse Bohr radius, the NLL correction 
vanishes. Furthermore, at v" = 0.86 GeV, the result becomes independent of v\ i.e., the 
NLL curve shows a local maximum corresponding to E^ s = 50 MeV, which we take as 
the central value of our estimate. The NLL curve also shows an impressive stability with 
respect to the hard matching scale variation in the physical range m c < < m^, as we 
observe in Fig. [2J The NLL curve has a local maximum at = 2.05 GeV, which we take for 
the numerical estimates. All this suggests a nice convergence of the logarithmic expansion 
despite the presence of the ultrasoft contribution where a s is normalized at the rather low 
scale v 2 jv h ~ 0.5 GeV. 

Let us discuss the accuracy of our result. For a first estimate of the error due to uncalcu- 
lated higher-order contributions, we take 6 MeV, the difference of the NLL and LL results at 
the local maxima. A different estimate can be obtained by varying the normalization scale 
in the physical range 0.8 < v < 1.4 GeV. In this case the difference with the maximum is 
14 MeV. Being conservative, we take this second number for our estimate of the perturba- 
tive error. Within the power counting assumed in this paper, the nonperturbative effects are 
beyond the accuracy of our computation and should be added to the errors. Following [9], 
we infer them using charmonium data. For an estimate we attribute the whole difference 
between perturbation theory and the experimental result, ~ 5 MeV, to nonperturbative 
effects. The nonperturbative contribution in heavy quarkonium is suppressed by the quark 
masses at least as l/(mim2m r ) and should be smaller for the charm-bottom bound state. 
We, however, take 10 MeV as a conservative estimate of the nonperturbative contribution 
to the HFS in B c . A further uncertainty is introduced by the error of a s (Mz). In Figs. [T] 
and [2] this is reflected by the yellow band, which is based on a s (M z ) = 0.118 ± 0.003. At 
the scale v" = 0.86 GeV it induces an uncertainty of MeV. 

To conclude, we have computed the HFS for a heavy quarkonium composed of quark 
and antiquark of different flavors in the NLL approximation by summing up the subleading 
logarithms a™ ln n_1 a s to all orders in the perturbative expansion. The use of the NRG 
stabilizes the result with respect to the v variation at the physical scale of the inverse Bohr 
radius and allow for solid first principle theoretical predictions. An explicit result for the 
two-loop soft anomalous dimension of the spin-flip potential is also presented. 

We predict the mass splitting of the vector and pseudoscalar charm-bottom mesons 

M(B*) - M(B C ) = 50 ± 17 (th) t\t ( 6a s) MeV (24) 

where the errors due to the high-order perturbative corrections and the nonperturbative 
effects are added up in quadrature in "th", whereas ll 5a s " stands for the uncertainty in 
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Figure 1: HFS for charm-bottom quarkonium as the function of the renormalization scale 
v in LO (dotted line), NLO (dashed line), LL (dot-dashed line), and NLL (solid line) ap- 
proximation for Uh = 2.05 GeV. For the NLL result the band reflects the errors due to 
a s (M z ) = 0.118 ±0.003. 



a s (Mz) = 0.118 ± 0.003. With improving statistics and precision of the B c data our result 
can be considered as a prediction for the B* meson mass. 
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Figure 2: HFS for charm-bottom quarkonium as the function of the hard matching scale 
v h in LO (dotted line), NLO (dashed line), LL (dot-dashed line), and NLL (solid line) 
approximation for v = 0.86 GeV. For the NLL result the band reflects the errors due to 
a s {M z ) = 0.118 ±0.003. 

Appendix 

The analytical results for the coeffcients fi and Ai of Eq. ( TH1) read (z = (a s (z/ p )/a s (z// l )) 1///3 °) 

f 3 = ^-4C A ) h = z 2/3 -3C A ? /g = ^-2C A ? /g = ^-2C A {q ^ _ ^ ^ 
f 7 = z 2 Po -C A ^ /g = ^o-(13C A )/6 ^ /g = Z p -2C A ? /iq = zPo+ C A ) 

/u = **\ /i2 = ^o /2 ( 2 -^' 1 + ^)' /is = **, / 14 = ^ln(2-^), 
/i 5 = z 3Ca , / 16 = ln(z) , f 17 = 1 , / 18 = In (2 - , (25) 
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+ (9C A - 8 ni T F )(5C A - 4mT F f(llC A - 4 ni T F )(C A + n{T F ) 

{-mAC A C%{C A + ni T F ) + 27C 2 A C 2 (7C A + 4n,T F ))(C A - 8n{T F ) 
8(5C A - 4n*T F )(C A - 2n{T F ){7C A - 2n{T F ){C A + n,T F )(9C A + 8n{T F ) [ M ' 



+ fi 



3C 2 



4(11C A - 4mT F ) 
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6C A (C A - 2C F ) 
(nC A -ArnT F ) ' 

19UC A C F (1SC A - 8mT F ) 27C A C F (3C A - An t T F 



13(5C A - A ni T F )(llC A - An{T F ) 2 {C A - 2n{T F ) (5C A - A ni T F )(llC A - An{F F ) 2 
27C A C 2 F {inC\ + mC A niT F - 416n 2 7f) 



+ /i. 
+ 



26(5C A - 4n i T»(llC A - A ni T F ) 2 (C A - 2n ; 7» 
-216C|(585C A - 554C A rnT F + 6AnfT 2 ) 



13(5C A - AmT F ){llC A - A ni T F ) 2 (C A - 2n{r F ) 
-hAC A C 2 F {32hC\ - 12Q8C A mT F + 2QAn 2 T F ) 



13(5C A - AmT F )(llC A - AmT F ) 2 (C A - 2n{F F ) 

2\W\C F + A32C 2 A C 2 F 1728C A C F + 8QAC 2 A C 2 F 

+ ji r 



(5C A - A ni T F )(llC A - An{T F ) 2 ^ (5C A - A ni T F )(llC A - A ni T F ) 2 

-86AC A C F (C A + nflp) + 27C 2 A C 2 F (7C A + A ni T F ) 
A(5C A - AmT F ) (C A - 2n t T F ) (C A + n t T F ) (9C A + 8n,7» [ M ' 

129QC\C F + A32C\C 2 (3C A - ni T F ) 
(9C A - 8n{T F ){hC A - AmT F ){llC A - An{T F ) 

-216C A C F (21C A - 8mT F ){5C A - 4n,7» - 1296C A C 2 F {AC A - 3n{T f 

+ ji r 



(9C A - 8 ni T F )(5C A - AmT F ) 2 (llC A - A ni T F ) 



F (,, A 2C A \ C F (C A + 2C F )(C A -8 ni T F ) 

' 2Fl I'' 4 " ft"' ~ l ) ^C A -An l T F ){9C A -An l T F ){2C A -n l T F ) { ° A + A ° F ^ 

_ / 2C A 2C A \QC A (C A -2C F ) 3C A (C A -2C F ) 
- B 1/2 2 - —, 1 + — - — — — — - + 



+ 



A) ' fa J (nc A -A ni T F ) i6(2C A - ni T F ) 

-3C A C F 



+ 



8(13C A - 8mT F )(19C A - 8n,7»(5CU - AmT F )(llC A - An{T F ) 2 {2C A - n x T F ) 
x (263641(7^ - 9\9\\AC A n{T F + 1071256C> 2 T£ - 556448C>fT| 

+ 131A56C A nfT F - 11264nfT|) 

27CUCJ. 



(9C A - 8n,7»2(19CU - 16n / T F )(5C A - 4n z 7» 2 (llC A - 4n z r F ) 2 



1 



x 7^ 9 T w 9r , (3644181C1 - 7690472C>7> + 3453968C> 2 T 2 

{(U A — 2nil F ){2u A — nil F ) v 

+3026560C>fT| - 3419648C>fT| + 1150976C A nf7| - 131072n?7|) 

3C 2 



+ 



16(19C A - lQn l T F )(9C A - 8n l T F ) 2 {l9C A - 8n,T»(5C A - An{T F ) 2 

X (11C A - An lTF ) 2 {7c!- 2n lTF )(2C A - n{r F ) ( 12488524839 ^ 

- 37966954860C>/Tp + 37940834480C> 2 T| - 13361 158A0C A nfT F 

- 27950404608C>fT^ + 25870953472C A nf T| - 1 14482053 12C 3 A rfT F 
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+2764505088C>[T£ - 343932928C A nfT| + 16777216nfT|) 

" -3Cf 

_(19C A - 16ra,7»(9CU - 8^Tp) 2 (5C A - 4^T» 2 (11C A - 4ra,7» : 



+ 



x 7^ 9 T L r ^-(62685009C A - 91230606C>,T F 

i'^A- ^niiF)\^A — ni±F) 

- 78455168C A n 2 T| + 2337725 YlC A n\T% - 176816384C>fT| 
+ 58415104C>fT| - 7979008C A nfT| + 2621UnjT F ) 

'F 



-3C? 



4(19C A - 16ra,7»(9CU - 8^T F ) 2 (5C A - 4n,T F ) 3 (llC A - 4n i 7» 2 

X (7C A -2n^(2C A -n^) ( 659490741 ^ " 1386410130^^ 

- 876382076C> 2 T| + 5528200720C> 3 T| - 74225 17824C>fT^ 
+ 5 15625 1904C A nf T F - 2102788096C A nfT| + 511131648C A n; 7 T^ 

- 69730304C A nfT| + 4194304n^T| 

-864C A C| - 432C A C F -3456C A Cf - 1728C A C| 



(5C A -4n,7»(llC i4 -4n,7»2 (5C A — 4niTp)(llC A — 4n;Tp) 2 



(26) 



with /i r = m r /{mi +m 2 ), B z (a,b) is the incomplete beta-function, and 2 Fi(a,b;c; z) is the 
hypergeometric function. 
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